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00 ABSTRACT 

We explore the microwave anisotropies at large angular scales produced by the emission from 
cold and large dust grains, expected to exist in the outer parts of the Solar System, using 
a simple toy model for this diffuse emission. Its amplitude is constrained in the Far-IR by 
the COBE data and is compatible with simulations found in the literature. We analyze the 
_C templates derived after subtracting our model from the WMAP ILC 7 yr maps and investigate 

on the cosmological implications of such a possible foreground. The anomalies related to the 
Q low quadrupole of the angular power spectrum, the two-point correlation function, the parity 

and the excess of signal found in the ecliptic plane are significantly alleviated. An impact of 
this foreground for some cosmological parameters characterizing the spectrum of primordial 

1 ^ i density perturbations, relevant for on-going and future CMB anisotropy experiments, is found. 

^sl Key words: Interplanetary Medium - (Cosmology): Cosmic Microwave Background - In- 

frared: Solar System - Submillimeter 
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O 1 INTRODUCTION AND TOY MODEL 

o 

I The thermal emission from Interplanetary Dust Particles (IDPs), called Zodiacal Light Emission (ZLE), is a well known source of foreground 
O in the (Far)IR sky. Recently |Maris et al.|p006] > discussed the use of the ZLE COBE model ( |Kelsall et al. |1998||Fix~sen & Dwek|2002| l below 



r- H 1 THz to derive predictions for surveys dedicated to Cosmic Microwave Background (CMB), such as WMAPNand Planck 
^ significant, contribution at frequencies t 857 GHz will be likely detectable by Planck. The simple scaling law of 



A weak, albeit 



Fixsen & Dwek 



i 2002 i 



predicts a small but not negligible contribution at 70-150 GHz, of particular interest in CMB studies. Theoretical estimates ( |Babich, blake| 
r> |& St einhardt 2007a b I suggest that the collective emission of Kuiper Belt Objects (KBOs) and other minor bodies could produce imprints in 
^ CMB surveys. 

CMB anisotropy maps derived from WMAP data are well consistent with the standard ACDM cosmological model. However, there 
are some intriguing deviations at large angular scales, such as those related to the angular power spectrum (APS) at low multipoles t's 
( |Bennett et al.|2003a||Hmshaw et al.|2003{|Spergel et al.|2003||Larson et al.|20T0)|Kim & Naselsky|2010a|b||Gnippuso et al.|2010[ >. Other 
anomalies regard the alignment of low multipoles (Tegmark, Oliveira-Costa & Hamilton 2003 , Gordon et al. 2005 , Copi et al. 2007 ; Huterer 
|2006[|Copi et al.|2007||Huterer|2006[|Gruppuso & Bur igana 2009 ; Gruppuso & Gorski 2010J and to the North-South asymmetry of the APS 
with respect to the ecliptic plane (Tegmark, Oliveira-Costa & Hamilton 2003; de Olive ira-Costa et al.|2004||Copi et al.|2"007l|Huterer|2 006). 
Moreover, Vielva et al. (20041 detected a localized non-Gaussian behavior in the southern hemisphere using a wavelet analysis technique, 
see also ( |Cruz et al.|2005) , and|piego et al. (2010) recently identified an intriguing excess of emission in WMAP 5 years ILC map along the 
ecliptic plane. Their results will be discussed later in more detail. 

If not produced by systematic effects (Tegmark, Oliveira-Costa & Hamilton 2003, Huterer 2006; Burigana et al. 2006, Gruppuso et al. 
2007) these features could be of cosmological origin (Gordon et al. 2005 , Ghosh et al. 2007) or effects produced by some foreground within 
the near Universe ( |Vale|2005||Cooray & Seto|2005||lnoue & Silk 2006a b, Rakic, Rasanen "&~Schwarz|2006") , the Galaxy ( |Frisch|2005) , or 
the Solar System ( |Schwarz etaT . 2004 , Copi et al. 2007). 
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17.5° 


2.905 


3.016 


3.362 


WM: 


17.5° 


1.211 


1.223 


1.256 


CM: 


35° 


2.698 


2.698 


3.100 


WM: 


35° 


1.194 


1.204 


1.234 


CM: 


70° 


2.064 


2.122 


2.303 


WM: 


70° 


1.135 


1.142 


1.162 



Table 1. Ratio between the quadrupole derived from the subtraction of the KBOE template from the WMAP 7 yr ILC map and the original WMAP map for 
the three WMAP Q, V, W bands and the considered models, using anafast. Models are the Cold Model (CM), with T = 30 K, r = 3 X 10~ 7 and the Warm 
Model (WM) with T = 60 K, r = 3 X 10~ 8 . 



It can not be excluded that some diffuse Solar System emission could contribute at low Ts at WMAP frequencies (Schw arz et al.|2004| 
|Copi et aT|2 007 ). Such foreground would exhibit planar symmetry with respect to the ecliptic (or some plane slightly tilted with respect to it). 
According to (Far)IR measures, the "classical" ZLE is the thermal emission of sub-mm IDPs, with orbital radii within 5 AU, characterized 
by a a; v 4 scaling below 1 THz I Fixs en & Dwek|2002 >. If larger particles exist, their emission would add to the classical ZLE, producing a 



signal at mm wavelengths, ranging between about 1/10 and 10 times the classical ZLE, with a less steep spectrum, as 1/v, as suggested by 
several models < |Backman, Dasgupta & Stencel|1995l|Stern|1996l|Yamamoto & Mukai|1998l|Moro-Martin & Malhotra|2002l[2003t . Such 
particles can be produced by the erosion of KBOs due to mutual collisions or by the erosion of incoming interstellar dust and they should 
be located beyond Jupiter. They would be then colder and possibly larger than the IDPs responsible for the classical ZLE. Impact detectors 
onboard interplanetary probes reveal that such particles indeed exist (Landgraf et al. 2002), but very little is known about their nature. 

Motivated by these considerations, we reconsider the problem of diffuse Solar System emission at cosmological frequencies, by in- 
troducing a toy model of emission from the KBO dust particles, that we will call in this paper KBO Emission (KBOE), and evaluating its 
implications for CMB observations. 

To explore the effect of the KBOE on CMB maps we considered a toy model in which the emission is confined within a region of width 
H kboe > symmetric with respect to the ecliptic plane, with constant brightness inside and zero outside. Owing to the planar and cylindrical 
symmetry of the model, all the coefficients a f m for the multipole expansion in ecliptic coordinates of the KBOE will be zero except for those 
with even £ and m = 0. Then KBOE (and likewise ZLE) will just affect the map components with even £. By denoting with I^ B0E (P), the 
KBOE for unit solid angle, integrated along a given pointing direction P, notable relations can be derived for the sky-averaged KBOE, /kboe, 
its variance over the sky and its quadrupole af B0E . Since in the model I™ 0E (P) is constant for P within an angle ±H KBOE /2 from the ecliptic, 
defining S = sin(//KBOE/2), we have /kboe = A KB oe<S , var(/™ 0E ) = A| B0E <S (1 - S) , and 

fl ™ OE = -AKBo E V5^S(l-<S 2 ) ; (1) 

here Akboe is the constant value of KBO dust emission within the region. The maximum variance is var(X) < ^kboe/4 and occurs for 
H kboe = 60deg, while the |«2,ol maximum is Iflo.ol < V20w/27Akboe ~ 1.53Akboe occurring for //kboe ~ 70.53 deg. Flux variations can be 
translated in term of brightness temperature variations with the usual conversion factor [<9B v (r)/<9r] 7 - =2 725K- 

Detailed information about the distribution of dust in the outer Solar System to fix //kboe an d A K boe is missing. A constraint on //kboe 
comes out from the distribution of inclinations of KBO orbits with respect to the ecliptic. KBOs are divided in two populations: classical 
KBOs with orbital inclinations within 10° about the ecliptic, and scattered KBOs with orbital inclinations within 40° about the ecliptic. In 
analogy with the distribution of dust from the erosion of Main Belt asteroids, we assumed that KBOs are a tracer of KBO dust and then of 
KBOE. Thus, we selected three possible representative values of //kboe = 17.5°, 35° and 70°. The region of the ILC map affected by the 
KBOE in these three cases is shown in Fig.[TJi (Note that this figure has been not included in the pubblished version for editorial reasons.). 

To constrain A K boe is a more complicated task, since it depends on details about grains such as their temperature, their size and their 
radial distribution, as well as their shape and their mineralogical composition. We rely on existing models predicting the possible sky- 
averaged KBOE brightness, leaving the development of a self-consistent model as a topic for a subsequent work. We choose as an example 
the model of Stern ( 1996), based on the production of dust in the KBO band by collisional disgregation of the KBOs. This model has the 
advantage to describe the sky-averaged KBOE in terms of two free parameters: the averaged dust temperature, T, ranging from about 10 K 
to 60 K and the sky-averaged optical depth, r < 10~ 6 . The model has been already constrained by the authors using IR data. Constrains 
come also from the microwave COBE/FIRAS data for the sky-averag ed ZLE, /f- E , {Fixsen & Dwek|2002| . Fig.fl]} compares the spectral 



energy distribution (SED) of the sky-averaged flux from KBOE, If- , from the Stern ( 1996) prescription with If . All those combinations 



of (T, t) values for which /J 0T is upper-bounded by the gray band in Fig.|TJj are acceptable. This is equivalent to take the (lf B0E /lf LE ) ratio 
at 350 GHz to be less than a 6.5. Fig. [I]) shows that models with dust temperature T = 10 K, 30 K and 60 K and r < 3 x 10~ 7 are allowed by 
COBE/FIRAS data, and they predict at 93.5 GHz /^ BOE in excess of up to ^ 80 times /f LE . At frequencies above 1 THz, KBOE is one to three 
order of magnitudes below the classical ZLE, and it does not display any seasonal effect, so that its contribution to the bulk ZLE emission 
seen in COBE/DIRBE is negligible. Note that in the case of the most extreme models allowed by the COBE/FIRAS data it is possible to 
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Figure 1. Panel a): the ILC map derived from WMAP 7 yr data. We displaied also with black lines the ecliptic plane and the three increasing regions where 
the KBOE is present in the case of our toy model, for three different increasing values of the height scale, //kboe = 17.5°, 35° and 70° Note that this figure has 
been not included in the pubblished version for editorial reasons. Panel b): Comparison of ZLE fluxes compatible with COBE/FIRAS data (Fixsen & Dwek 
2002) and a set of possible models of KBOE (Stern 1996). The black solid line shows the ZLE derived from the best fit model to COBE/FIRAS data and 
extrapolated to lower frequencies. The gray band represents a sketch of the allowed region obtained from the error bars in Fixsen & Dwek (2002). The blue 
dashed lines display four different models of KBOE corresponding to different values of r and T. The resulting fluxes, sum of KBOE and ZLE, are represented 
by the red solid lines. Note that the classical ZLE (estimated on the basis of COBE data) is negligible in practice at WMAP frequencies, whereas KBOE might 
not be ignored. 



predict a contribution to the quadrupole (see Eq. |T]l) of up to ss 10 2 //K 2 , to be compared with the expected ZLE contribution less than 



0.02 [iK 2 and the CMB quadrupole moment AT"cmb/=2 ~ 200 yuK 2 ( |WMAP Explanatory Supplement| ( |2006) , see also footnote 1). 

How massive must be the debris disk to produce the KBOE? An estimate of the required Kuiper Belt mass to be converted in dust, Atj U8 t, 
can be readily derived from t. From dynamical models of dust propagation and the distribution of objects in the Kuiper Belt it is possible 
to infer that Kuiper Belt Dust should be confined within R in « 30 AU and R oul « 50 AU. Since T oc 1 / Vfl the grains could be considered 
approximately isothermal. As a further approximation it is possible to assume the numerical density of the dust n to be uniform between R m 
and R out . Assuming grains as spheres, with a typical diameter, a, unit emissivity efficiency, scattering cross-section given by their geometrical 
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Figure 2. APS at low multipoles as derived by the WMAP team analyzing the ILC 7 yr map (blue solid line) with its lcr errors (green solid lines); APS of 
the ILC 7 yr map derived using Anafast (black solid line); APS of the ILC 7 yr map after the subtraction of our KBOE template with //kboe = 70°, 35° and 
17.5° (black dotted-dashed line, dashed line, and dotted line, respectively). The dust model with r = 3 X 10~ 7 and T = 30 K (CM) and the V band is here 
considered. The red solid line shows APS of the best fit ACDM model for the WMAP 7 yr map. 



cross section, then n « 4r/[na 2 (R out - R in )] and Mi mt ~ (87r/9)rpa sin(//KBOE/2)[^ 3 ut - R 3 m ]/[R ut ~ ^inL Fragments should be composed 
largely of silicates, ices, and carbon compounds with densities p ranging between 1 and 3 gr/cm 3 . For t < 3 x 10~ 7 , //kboe :S 70°, and 
a«l cm, /Hdust varies between 6 x 10 23 gr and 18 x 10 23 gr, i.e. less than 0.14 Pluto masses. 



2 ANALYSIS 

The toy model outlined in Sect. [T] is exploited to assess its possible implications for CMB anisotropy statistical estimators and thus for 
cosmological models and parameters, under the assumption that such a kind of component present in the microwave maps, but it is hidden in 
the currently available CMB component maps. 

Template production - According to the above assumption, we subtract our toy models from the available CMB anisotropy map. We will 
present here results obtained exploiting some KBOE toy model templates characterized by //kboe = 17.5°, 35°, and 70° for two combinations 
of averaged dust temperature and sky-averaged optical depth: r = 3 x 10~ 7 , T = 30 K (cold model, CM) and r = 3 x 10~ 8 , T = 60 K (warm 
model, WM). These templates are computed in the ecliptic frame at the centers of the WMAP Q, V, and W frequency bands and generated 
at a resolution defined by the HEALPix (Gorsky et al. 2005 1 parameter N side = 256 (i.e. with a pixel size of about 13.7'), appropriate 



to our large scale analysis. They are then convolved with a Gaussian symmetric beam with FWHM of 1° to match the beam resolution of 
the ILC 7 yr map released by the WMAP team and then transformed into Galactic frame. We have then subtracted the proper monopole to 
each template. At each frequency, we subtract our KBOE templates from the WMAP ILC 7 yr map (previously degraded at N sjde = 256) to 
produce new CMB anisotropy maps cleaned by the adopted KBOE contribution. 

Angular power spectrum - We extracted the whole-sky APS of the ILC 7 yr map released by the WMAP team and of the maps cleaned 
by the KBOE contribution described above using the publicly available HEALPix routine anafast, applying 3 iterations and subtracting 
monopole and dipole. The APS derived for ILC 7 yr map is compared with that released by the WMAP team. The results are shown in 
Table [T] and Fig. [2] for the low multipoles of interest here. As evident, no significant difference is found between the APS derived with 
anafast and with the finer approach applied by the WMAP team, thus probing that the results simply derived with anafast are accurate 
enough for the aims of this work. Table[T|shows that the subtraction of the KBOE template from the ILC map in the case of the CM implies 
a larger enhancement of the CMB quadrupole, while the effect, although present, is not so remarkable in the case of the WM. Also, this 
enhancement increases with the decrease of Hkboe- These results do not depend significantly on the considered frequency channel because 
of the considered KBOE spectral shape. Thus, we will focus in particular on the V band, where Galactic foregrounds are minimum at large 
scales, and on the CM. In order to evaluate the probability that CMB quadrupole phases anti-correlate with the "deterministic" ones of 
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Figure 3. Two-points correlation function computed for maps at HEALPix resolution N$nie = 16. The solid line displays the average of 10 5 MC realizations 
of CMB anisotropy maps extracted from the WMAP 7 yr best fit ACDM model APS shown in Fig. [2] The red dotted lines show the corresponding lcr level 
fluctuations of the MC simulation. The blue dashed line refers to the ILC 7 yr map. The green dotted-dashed line refers to the map derived subtracting from 
the ILC 7 yr map one of our KBOE model, namely that obtained for the CM with //kboe = 17.5°. 



KBOE, we have performed a simple Monte-Carlo simulations. We extracted 10 4 random quadrupole realizations from the best fit of ACDM 
model obtained by WMAP 7 yr data and, for each of them, we have added our template and evaluated the modified quadrupole. We found 
that for about one out of three of realizations the band power, £(£ + l)C{/(2n), of the quadrupole decreases from its primordial value, and, for 
example, for a 11.5 — 13.4% of the events it is lowered more than 300 fiK 1 depending on the specific CM considered^] These numbers do 
not depend strongly on the considered value of /Ykboe- This shows, reversely, that the probability of observing an ILC map whose quadrupole 
turn to be increased, once the KBOE template is removed, is not negligible. Fig. [2] reports the APS up to £ = 20. Note that only the power 
of even multipoles is significantly affected, as expected since the geometrical symmetry of the KBOE component. For //kboe = 70° only 
the quadrupole turns to be significantly amplified, while for #kboe = 35° the power enhancement is important also at I = 4 and 6 and not 
negligible at t = 8. Remarkably, for Hkboe = 17.5° the power increase is important for the even multipoles up to I = 16 (see also the 
following discussion on Parity). Note that the overall effect, when KBOE template is subtracted from the ILC map, is an increase of the 
resulting CMB power at even multipoles. This is for geometrical reasons, but related to the relative phases, or orientations, of a Cm of ILC and 
of KBOE patterns. 

Correlation function - The two point correlation function is defined as C{6) = (T(hi)T(h2))r, v rn=cos(e) . where the symbol {..■)h 1 h 2 =cos(S) stands 
for the average over an ensemble of realizations in which the directions hi and h 2 form an angle given by 6. In Copi et al. ( 2007 ) it is shown 
that at angular scales greater than about 60° the WMAP data occur in 0.025% of realizations of the concordance model (see also Co pi et ah] 
1 2007^ ; |Hajian||2007| ); |Copi et al.|j2010) and references therein). In Fig.[3]we show C{9) (black line) obtained through 10 5 MonteCarlo (MC) 
random realizations extracted from the WMAP 7 yr best fit ACDM model whose APS is displayed in Fig. [2] These random realizations 
properly take into account the anisotropic white noise level of the V band channel of WMAP 7 yr data. The resolution we have considered 
for this analysis is given by the HEALPix parameter N^e = 16 and the angle 8 has been binned with a size of 5°. The red dotted lines show 
the corresponding lcr level fluctuation of the MC simulation. The blue dashed line refers to the ILC 7 yr map. The green dotted-dashed line 
refers to the map derived subtracting from the ILC 7 yr map one of our KBOE model, namely that obtained for the CM with H KB oe = 17.5°. 
It is interesting to note how well the green line is (almost) always well consistent within lcr level with the ACDM model whereas the blue 
line (i.e. the WMAP 7 yr ILC map) is often out of the lcr contour, especially at large angular scale. This improvement is typical of our KBOE 
templates for the CM, while the templates for the WM do not change much the WMAP 7 yr ILC map correlation function. 
Cosmological parameters - We quantify here the impact of the removal of the KBOE emission on cosmological parameters. The increase of 
the amplitude at low multipoles has interesting consequences on some of the cosmological parameters. For this purpose, we use the WMAP 
7 yr likelihood code publicly available with the option of a pixel likelihood code at low resolution and substitute the WMAP 7 yr ILC map 
with the KBOE cleaned one. The Markov Chain Monte Carlo package CosmoMC (Lewis & Bridle 2002 1 is connected with this modified 
likelihood code. As shown in Fig.|4]the removal of a Solar System contamination (namely the CM with H KB0E = 35°) modifies slightly 



Considering the fiducial ACDM model as a reference, the decreasing of at least 300 jiK 1 occurs for ss 30 - 36% of the events. 
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Figure 4. Two dimensional marginalised probability distributions for cosmological parameters by removing (black lines) or not removing (red lines) the 
KBOE template for the CM with //kboe = 35° (curves are the 68% and 95% confidence level). To the left (right) the plot for A s (n mn ) vs n s for a standard 
ACDM model (including running of the scalar spectral index). 



the amplitude and the spectral index of scalar perturbations (by choosing a pivot scale k, = 0.002Mpc ) in a standard ACDM scenario 
characterized by 6 cosmological parameters (plus the amplitude of the residual Sunyaev-Zeldovich effect): as expected by Fig. [2] a slightly 
larger value for the amplitude of primordial fluctuations is preferred by the larger amplitude of the temperature anisotropies at low £'$. We 
find a little impact of the removal of this foreground on the constraints on the spatial flatness of our Universe. A larger effect is expected in 
extensions of the ACDM model in which the small value of the quadrupole plays a relevant role: Fig. [4] shows how smaller values of the 
running of the scalar spectral index n mn = dn j(d\nk) are preferred once the KBOE is removed (in this last case we have adopted an optimal 
pivot scale k, = 0.017Mpc~'). It is important to note how the change in the mean value of « mn we obtain is 0.023 after removing the KBOE 
template. Such difference is comparable with the Planck sensitivity to M m n: it is therefore important to assess foreground contamination on 
large scales to infer on the primordial spectrum of density perturbations. Also, since correcting for KBOE the low multipole tail in TT is not 
so suppressed, including tensor perturbations in the analysis, we find that, the upper limit on tensor-to-scalar ratio of primordial perturbations, 
r TS , is relaxed to r TS < 0.42 (at 95%CL), compared 0.36 given in |Larson et al.| ([2010). Note that this effect is due to the larger relevance of 
temperature information with respect to polarization in WMAP data. 

Parity - One of the large scale anomalies is the Parity anomaly. For all-sky maps it is customary to expand the CMB temperature fluctuations 
T(h) in terms of Spherical Harmonics Y( m (h) where h is a direction in the sky, namely depending on a couple of angles (0,(f>), T(h) = 
YiCm a im Yi m (n) = Yit Te («) , where in the last equation we have implicitly defined the maps T[(h) for each angular scale £. Since under 
reflection (or Parity) symmetry n — » —h the Spherical Harmonics behave as Y (m (n) — > (-l) f Y Cm {n) , and the map T e (n) has even parity for 
even t and odd parity for odd I. 

At large scales, we expect a Sachs-Wolfe plateau in terms of TT APS, therefore we roughly forecast the same amount of power in even 
and odd maps. It has been recently proved that for WMAP data at the large scales ( Kim & Naselsky 2010a b I, precisely in the range of 

6 = 2— 22, the total power coming from even £'s is unlikely smaller than the total power present in the odd £'s. The probability related to 
this event is as low as 0.4% for the WMAP 7 yr data Kim & Naselsky (2010b) (see Gruppuso et al. ( 2010) for a temperature and polarization 
joint analysis). The analysis of this power asymmetry was first proposed in |Land & Magueijo (2005b]( to test the presence of foreground 
residuals since templates of dust, free-free or synchrotron emission possess a large Parity asymmetry. Also the KBOE template we propose 
in this paper, as well as the standard ZLE, is highly Parity-asymmetric. Therefore it is worth to study its possible impact on the WMAP ILC 

7 yr map, of course under the assumption that a residual of the kind considered in the current paper, is present in the ILC map. 

We consider the same estimator of Kim & Naselsky 1 2010a b i defined as r = P + /P~, where = £(£ +1) Ccl2n ; here with 
Yic+i- we mean the sum over even or odd £'s respectively in the considered range (here 2 - 22^J. We have extracted 10 5 random maps from 
the best fit model of WMAP 7 yr working at the HEALPix resolution N si{te = 64 (i.e. all-sky maps of 49152 pixels). We have computed the 
r estimator for each random map and we have built the probability distribution function (pdf) for r, shown in Fig. [5] Note that the pdf for r 
does not peak around the value 1, but slightly larger since for the chosen I -range, there are more terms at numerator than in the denominator. 
Vertical lines in Fig. [5] represent the values of r for the considered maps. Black vertical line (on left) stands for the WMAP 7 yr ILC map. The 



We consider the same range of |Kim & Naselsky] j2010at where it was shown the maximum of the anomaly. See also |Bennett et al. ]l [2010fr for comments 
about "a posteriori" selections. 
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Figure 5. Parity anomaly of the estimator r = P + /P_ as defined in the text with f max = 22. The histogram (in red) displays the distribution of r obtained from 
10 5 MC realizations. Vertical lines correspond to the maps considered in this work: the black solid line (on the left) refers to the ILC 7 yr map; colored solid 
lines refer to the CM; colored dashed lines refer to the WM. Green, blue, and yellow lines are for //kboe = 17.5°, 35°, and 70°, respectively. 



probability to have this value is as low as fj.91%0 When we remove our KBOE template from the maps we have the colored vertical lines. 
Green, blue, and yellow lines stand for H KB0E = 17.5°, 35°, and 70°, respectively. Colored solid lines refer to the CM, colored dashed lines 
refer to the WM. Note that the probability associated to the green line (i.e. CM and H KB oe = 17.5°) is 77.1% and the probability associated to 
the blue line (i.e. CM and H KB oe = 35°) is 10.6% . The Parity anomaly is removed when these templates are properly taken into account. In 
the definition of the Parity estimator in the current paper we adopted £ max = 22, just for simplicity and since for that ( max value the anomaly is 
remarkable. However, as shown in Kim & Naselsky (2010a b I; Gruppuso et al. ( 2010), in the WMAP TT spectrum there is a whole multipole 
range, rather than a single l max value, where the WMAP 7 parity anomaly holds. This dims significantly the case for posterior biasing. In 
analogy to what performed for the quadrupole, we have implemented a simple Monte-Carlo simulation in order to evaluate the probability 
that CMB even multipoles phases (up to l max = 22) anti-correlate with the "deterministic" ones of KBOE. We extracted 10 4 random CMB 
realizations from the best fit of ACDM model obtained by WMAP 7 yr data and, for each of them, we have added our template and evaluated 
the modified spectrum. Hence, considering only the CM, we have found that the probability of lowering the estimator r is 24.7%, 7.05% and 
0.07% for -Wkboe = 70°, 35° and 17.5°, respectively. Similarly, we computed the probability of lowering the "even powers" up to £ max = 22 
for some fixed threshold. For instance, for Hkboe = 70° the probability to decrease the total even band powers of at least 300/j.K 2 with respect 
to the realization (or to the fiducial ACDM model) is 7.46% (43.46%), whereas it is 2.67% (23.26%) for H KBO e = 35° and 0.01% (0.67%) 
for //kboe = 17.5°. This shows that for //kboe larger than few tens of degrees the anti-correlation probability turns out to be not negligible. 
These numbers show that, contrarily to the quadrupole case discussed above, in this case the dependence on H KBOE is important. 
Alignments - We have seen how the removal of the KBOE template (for the CM) can reduce the unlikeliness of some estimators, as the 
C( and those related to them. Our preliminary tests carried out using a public code for the multipole vectors decomposition p] by |Copi et al. 



( |2004} show that this does not happen for the estimators based on ac„, phases by means of which alignment anomalies are detected. 
Ecliptic excess in WMAP data - In |Diego et al"j ( |2010fr , the authors look at the linear combination V + W - 2Q of WMAP bands. In 
their analysis, the WMAP bands have been cleaned of Galactic contamination and the bright point sources have been masked. The linear 
combination V + W - 2Q gets rid off the CMB signal completely and should contain a linear combination of the instrumental noise in 
each band, beam effects and residual Galactic and extragalactic foregrounds not accounted for in the cleaning (and masking) process. After 
masking out the Galactic plane, the authors found a significant large scale signal aligned with the ecliptic plane with a grey-body spectrum 
(see Fig. 1 in Diego et al. (2010)). A Galactic origin for this signal was ruled out due to the high latitudes where it was found. An extragalactic 
origin was also ruled out due the highly anisotropic distribution of this signal. In that paper it was suggested a possible origin in the ZLE 
but extrapolations from actual measurements from COBE/FIRAS predicted a much smaller signal than needed to account for the excess in 
the ecliptic plane. The present paper offers a fresh alternative to explain the excess by assuming a different population of larger dust grains. 
We have combined the signal expected from our KBOE models with the instrumental noise and beam of WMAP and built templates of the 
V + W - 2Q maps. In Fig.|6]we show the predicted signal from the three CM models. The comparison with the signal detected by Dieg o et af] 



6 This value is slightly larger than what quoted in |Kim & Naselsky|j2010a) that is 0.86%. This might be due to the details of the MC, like number of 
simulations, resolutions of maps and the considered FWHM. 

7 See also http://www.phys.cwru.edu/projects/mpvectors/ 
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Figure 6. Predicted signal in the combination V + W - 2Q from the CM with Hkboe = 70° (bottom), Hkboe = 35° (middle), and Hkeoe = 17° (top). The 
map units are /.(K and it has been smoothed with a 7° Gaussian. Compare this plot with Fig. 1 of Diego et al. (2010) based on WMAP data: the model with 
Hkboe = 70° seems to agree better with the observations. 

l |2010fr (see their Fig. 1) shows that the model with Hkboe = 70° seems to agree better with the observations than the others. The areas near 
the Galactic plane show no evidence of excess appearing in WMAP data, that could be due to residual synchrotron emission. Note that while 
KBOE should be positive in V + W - 2Q , the synchrotron signal shows up as a negative signal in the combination V + W — 2Q. 

3 CONCLUSION 

We exploited a simple toy model for diffuse emission from cold and large dust grains expected to exist in the outer part of the Solar 
System and requiring only a modest amount of mass, with signal amplitudes, constrained in the Far-IR by COBE data, compatible with 
simulations existing in literature. We have produced and analyzed templates derived subtracting our toy model from WMAP ILC 7 yr maps 
to investigate on cosmological implications of such a foreground. We find that the anomalies related to the low quadrupole of the angular 
power spectrum, the two point correlation function, and the Parity are significantly alleviated. No significant impact has been found on low 
multipole alignments; of course, this does not exclude that more detailed models constructed in this framework could have consequences for 
these estimators. 
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It is interesting to note that there might be a relationship between the lack of effect of our model on the alignment of the low multipoles 
and the low amplitude of the quadrupole measured by WMAP. If the CMB intrinsic quadrupole is anti-correlated (i.e. with a similar pattern 
but opposite sign) with the ecliptic plane, any signal originating in our Solar System (like those considered in this paper) would compete 
with the CMB quadrupole, reducing its measured amplitude. On the other hand, the subtraction of the Solar System signal to the measured 
quadrupole would not change the orientation of the resulting signal since the subtraction operation transforms the anti-correlated signal in 
a correlated one and both patterns add up. A visual example of this can be found for instance in Fig. 9 of Diego et al. (2010) where also 
a hypothetical signal that traces the ecliptic plane is assumed. The same authors also demonstrate how signals in the ecliptic plane affect 
(significantly) only the even multipoles. Consequently, the octupole (I = 3) would remain basically unchanged. The KBOE is able to explain 
the excess of signal found in Diego et al. ( 2010), in particular for cold dust particles and for a high H KB oe- 

Finally, we showed that this foreground has an impact for some cosmological parameters characterizing the spectrum of primordial 
density perturbations, relevant for on-going and future CMB anisotropy experiments. 

Clearly, the model needs to be improved in the future and tested against the data from the Planck satellite. Our analysis shows that it 
will be relevant to include, or at least to test, such component in the analysis of microwave anisotropy data at large scales. 
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